ON THE NUMBER OF ZEROS OF MELNIKOV FUNCTIONS 



SERGEY BENDITKIS, DMITRY NOVIKOV 

Abstract. We provide an effective uniform upper bond for the number of 
C ) ' zeros of the first non-vanishing Melnikov function of a polynomial perturba- 

Cn ' tions of a planar polynomial Hamiltonian vector field. The bound depends on 

degrees of the field and of the perturbation, and on the order k of the Melnikov 
^J i function. The generic case k = 1 was considered by Binyamini, Novikov and 

Yakovenko ( BNYIO]). The bound follows from an effective construction of 

the Gauss-Manin connection for iterated integrals. 
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1. Introduction 



1.1. Infinitesimal Hilbert 16th problem. The second part of 16th Hilbcrt prob- 
lem asks: How many limit cycles may have a planar polynomial vector field? The 
question has a long history, and was at the origin of several theories, see [102 ). The 
most remarkable achievement, Ecalle-Ilyashenko theorem, claims that the number 
of limit cycles is finite for any individual vector field, see I,E92;. However, exis- 
tence of a uniform upper bound for this number even for quadratic vector fields is 
an open problem. 
^ .— ' A weaker form of the same question concerns perturbations of Hamiltonian vector 

f-^ . fields. Let H(x,y) be a bivariate polynomial (further called Hamiltonian). The 

corresponding Hamiltonian system can be written in Pfaffian form as 

O dH = 0. (1.1) 



Consider its polynomial perturbation 

dH + euj = 0, where w = P(x,y)dx + Q(x,y)dy, P, QeR[i, y], (1.2) 
ande 6 (tt 1 ^). 



H \ Consider a nest of cycles {S t C {H = t},t £ [a,b] C R} of ([LIT) . We ask how 

many limit cycles of (|1.2p converge to this nest as e — > 0. 

It is easy to see that closed trajectories St that survive after the perturbation 
should produce zero value of the Poincare integral (aka first Melnikov function) 



/ = I(6t,u) = <t u, 
J s t 

the so called Poincare- Andronov-Pontryagin criterion, see |IY1 §26A]. Therefore 
estimates on the number of zeros of this so-called Abclian integral have direct 
relation to the Hilbert 16th problem. Binyamini, Novikov, Yakovenko studied the 
case of non- conservative perturbations, namely, when the Poincare integral does 
not vanishes identically. 
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Theorem 1. 'BNY101 Assume that 1^0 for the nest of cycles of 1 1.1)) . Assume 
that degw < degiJ. Then the number of cycles St providing the zero value of 
Poincare integral is at most 2 2 , where P(n) is some explicit polynomial of 

degree at most 60. 

This upper bound serves also as an upper bound for the cyclicity of an open 
nest of the limit cycles (which is defined as a supremum of cyclicities of all closed 
subnests of the open nest, see e.g. |GN10j ). 

For generic Hamiltonians identical vanishing of / implies exactness of w (again, 
assuming degw < dcg-ff), so the perturbation remains integrable, see [169) . How- 
ever, for degenerate Hamiltonians one has to consider Melnikov functions of higher 
order. 

1.2. Melnikov functions and the main theorem. 

Definition 1. For a cycle S of (|1.1|) choose a transversal a with coordinate z chosen 
in such a way that S intersect a at z = 0. 

Denote by A : a — > a the holonomy map of cycle 7 considered as a function of 
the parameters h, e. Being analytic function of its arguments, A can be expanded 
in the converging series 

A(z, e) = z + eA'hiz) + ■■■ + e K M K {z) + ..., (1.3) 

where Mk{z) are real analytic functions defined in some common neighborhood of 
the origin z = 0. The function Mj~ is called fc-th Melnikov function. 

Assume that the first nonzero function M^{z) has N isolated zeroes (counted 
with their multiplicities) in the closed interval {\z\ < p}. 

Proposition 1.1. [1Y , Proposition 26.1] There exists a small positive value r > 
such that for all \s\ < r the foliation (|1.2[) has no more than n limit cycles 
intersecting a at {\h\ < p\. 

Our main result provides an upper bound for the number of isolated zeros of the 
first non-zero Melnikov function. 

Theorem 2. The number of isolated zeroes of the first non-zero Melnikov function 

Mk is bounded by 2 , where n + 1 = degH, d — degw, and the absolute 

constants in 0(1),0(K) can be explicitly computed. 

This bound is certainly not exact, and construction of lower bounds is a difficult 
problem, unsolved even for Abelian integrals. 

Note that the order K of the first non-zero Melnikov function cannot be easily 
bounded in terms of degree of H: this problem includes, as a particular case, the 
center-focus problem. 

1.3. Iterated integrals and algebraic motivation. It is well-known, see [G05I 
IIYj , that Mk can be represented as a linear combination of so-called iterated inte- 
grals of order at most K. 

Definition 2. Let 7(s) : [0, 1] — >• C 2 be parameterization of a curve 7 C C 2 . For a 
fc-tuple of forms Wi,...,Wfc £ A : (C 2 ) we define the iterated integral as 

7 <« = / (I 1 (■■■(/ V«*) 7^-r) ...) TV- 
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Iterated integrals were extensively studied from various points of view, see e.g. 
|Ch|IHPMN08] . Our goal is to investigate their oscillation properties. Let us choose a 
straight line as a transversal to the nest of cycles. Iterated integrals define functions 
on this transversal: to any point p of the transversal corresponds the value of the 
iterated integral over the cycle of the foliation passing through it, with p being the 
initial point of the path of integration (note that, unlike the Melnikov function, the 
iterated integrals do depend on the choice of the initial point of the cycle) . 

The main step of the proof of Theorem [5] is an explicit construction of a mero- 
morphic flat connection whose horizontal sections are given by basic iterated in- 
tegrals (see (|3.6|) for definition), a higher order analogue of the Gauss-Manin con- 
nection for Abelian integrals. We prove in Section |4] that this connection belongs 
to the class of connections considered in the paper of of Binyamini, Novikov and 
Yakovenko [BNYlOj . see the next section for formulation of the result. Estimates 
on the complexity of the connection, proved in Section [31 allow to apply their main 
result not only to linear combinations of basic iterated integrals, but also to their 
combinations with coefficients polynomially dependent on z from (|1.3j) . In Section[5] 
we represent Mk in this form. 

2. NON-OSCILLATION OF HORIZONTAL SECTIONS OF MEROMORPHIC 

CONNECTIONS 

In this section we briefly recall the main result of |BNY10j . Let 17 be a rational 
I x Z-matrix of rational differential 1-forms on a complex manifold M , with a singular 
locus S. It defines a connection 

dX = Q ■ X (2.1) 

on trivial vector bundle M x C l . We denote by S the singular locus of the connec- 
tion. 

2.1. Regular integrable connections. 

Definition 3. The form fi is integrable or locally flat if dfl — O A O = 0. 

This condition is equivalent to local existence of a basis of horizontal sections of 
(|2.1[) near each nonsingular point a ^ E. 

Definition 4. The Picard-Fuchs system (12.11) (and the corresponding matrix 1- 
form Q) is called regular at a G M, if for any germ of a holomorpic curve 7 : 
(C,0) — > (M,a) the pull-back of the connection to (C,0) has a regular singularity 
at the origin: 

VOO 3p = P (C)eR ||X( 7 ( s ))|| ±1 = 0(| s rP) (2.2) 

as s — > in the sector {args| < C}. 

Connection is called regular on M if it is regular at each point a G M. 

Regular connections remain regular after pull-backs, push- forwards, (semi)direct 
products etc., see (D) . 

2.2. Quasiunipotent connections. 

Definition 5. For a point a € M a small loop around a is a closed path 7, such 
that exists a mapping {\z\ < 1} — > M which maps to a, {|z| = 1} to 7 and such 
that the image of {\z\ < 1} \ {0} is disjoint with S. 
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Recall that an operator is called quasiunipotent if all its eigenvalues are roots of 
unity, i.e. belong to exp(27rzQ). 

Definition 6. The integrable form Q is called quasiunipotent at a point a £ M, 
if monodromy operator associated with any small loop around a is quasiunipotent. 
The system is (globally) quasiunipotent, if it is quasiunipotent at every point of 

CP". 

In general, this does not mean that every monodromy operator associated to fi 
is quasiunipotent. 

Theorem 3. (Kashiwara theorem JK81]). A regular integrable system that is qua- 
siunipotent at each point outside an algebraic subset of codimension 2, is globally 
quasiunipotent. 

2.3. Degree of rational function. We define degree of a rational function to 
be the minimum of sums of degrees of numerator and denominator over all its 
representations as a ratio of two polynomials. Degree of the form is defined in such 
a way that the operator d has degree 0. 

2.4. Notion of size. In this work, similar to |BNY10j . we are studying various 
objects, like matrices, functions, differential forms, defined over Q - field of rational 
numbers. To obtain quantitative characteristics of these objects, we need to use 
the notion of size, or complexity of the objects. 

Definition 7. The norm of a multivariate polynomial P £ C[zi, ■ ■ ■ , z n ], P(z) = 
Tli a c a z a (in the standard multiindex notation) is the sum of absolute values of its 
coefficients, |P|| = ^ a \c a \. Clearly, this norm is multiplicative, 

\\PQW < \\p\\ ■ IIQII 

Definition 8. The size S(P) of an integer polynomial P 6 Z[zi, . . . ,z n ] is set to 
be equal to its norm, S(P) = ||P||. 

The size of a rational fraction R £ Q(^i, • ■ • , z n ) is 

S(P) = min{||P|| + HQII : R = P/Q: P,Q G Z[z u . . . , z n ]} 

The size of a (polynomial or rational) 1-form on P m or on P m x P 1 defined over 
Q, is the sum of sizes of its coefficients in the standard affine chart C m . 

The size of a vector or matrix rational function (resp., 1-form) defined over Q, 
is the sum of the sizes of its components. 

Note that, unlike polynomials, for rational functions we have only 

deg(^Pi) <2^degP, 

1 (2.3) 

R l ) <(n + l)IJs(i2i). 



2.4.1. Counting number of zeroes of the solution. Let Q be a rational / x Z-matrix 
1-form of degree d on the product CP m , and consider the restriction of the cor- 
responding Picard-Fuchs system (|2.ip to some line I = CP 1 C CP m . We are 
interested in the number of zeros of a linear combinations of entries of the funda- 
mental matrix of (|2.ip . In general, restriction of the fundamental matrix to this 
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line produces a multivalued matrix function on £ \ E, so to count zeros one should 
choose a simply connected domain in £ \ E. One can easily see that the geometric 
complexity of the domain should be taken into account. 

Definition 9. We denote by J\f(£) — Af(£l\e) the supremum over all constant 
matrices B and all triangles T lying in £ \ E of the number of isolated zeroes of the 
function Tr BX in T. 

Theorem 4. 'BNY10, Theorems 7,8] Let U be a rational I x l-matrix 1-form of 
degree d on the product <CP m x CP 1 . Assume that Q is integrable, regular and 
quasiunipotent, and its size is s = S(fl). Then 

W S CP 1 c CP m x CP 1 N{£) < s 2 C(d ' 4m)5 
for some universal constant C . 
3. Construction of Gauss-Manin connection for iterated integrals 

3.1. Base spaces: notations. Let C n +i[x,y] denote the space of all bivariate 
polynomials of degree at most n+1. We will denote the points of its projectivisation 
PC n +i[x,y] by A. In standard coordinates H — Eo<i+j<n+i ^ij^V 3 - By A we 
denote the tuple of all elements of A except the last one, Aoo- 

An important role plays the space C n +i[x, y] of polynomials vanishing at the 
origin, of dimension smaller by 1. The tuples A parameterize the points of its 
projectivisation PC n +i[x, y]. 

3.2. Gauss-Manin connection for Abelian integrals. 

Definition 10. Let H £ C[x,y] be a polynomial of degree n+1. Petrov module 
Ph is a C[i]-module defined as quotient space 

A 1 
H ~ dH ■ A + dA° 
of polynomial 1-forms over a space of relatively exact forms / • dH + dg, where /, g 
are polynomials. 

Proposition 3.1. 'IY, Theorem 26.21] The set of all Morse-plus polynomials H 
for which the forms u)ij — x z y 3 dx, 1 < i, j ' < n form a basis of Ph over C[h] is 
a Zarisky open subset in PC n +\[x,y]. 

Forms Wy provide a convenient trivialization of homological Milnor bundle over 
PCn+i[i,j/]. The Gauss-Manin connection in this trivialization can be written ex- 
plicitly and links the main result of [BNY10 to Infinitesimal Hilbert 16th problem. 
Let us formulate this result. 

Let H be a polynomial satisfying conditions of Proposition 13.11 such that the 
affine curve r# = {H — 0} C C 2 is smooth. Choose a point po £ Tjj- ^h 
is a Riemann surface of genus " 2 ~ with n+1 removed points. Therefore its 
fundamental group t^i{YhtPo) is a free group in N = n 2 generators. 

Choose Si, ...Sn £ F in such a way that their homology classes form a basis in 
Hi(T H ,Z). 

Theorem 5 ( |BNY10j ). The matrix 

Si = if ui\ , (3.1) 



<5fc J k,l=l 
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where u>i is an enumeration of the set of basic forms B = {u>ij = x l ~ x yi dx,i,j — 
1, ...,n}, is non- degenerate. Moreover, Si = Si(H) is the matrix of fundamental 
solutions of the Picard-Fuchs equation 

dSi = fiiSi, (3.2) 

which is defined over Q and has the size s — S (O), dimension £ and the degree 
d = dcg Q explicitly bounded from above as 

s<2 Poly (™\ d<0{n 2 ), £ = n 2 . (3.3) 

Using these estimates and Theorem [4] one get the main result of [BNYlOj . 
Our goal is to generalize this construction for iterated integrals of length K > 1. 
To this end we will need more detailed results. 

Proposition 3.2. Let 9 be a polynomial one-form of degree d on C 2 , and assume 
that it it is defined over Q(A), is of degree v in A and of size s. Denote by X — 
A \ {Aoo} the tuple of all coefficients of H except the first one. Then one can write 
a decomposition 

N 

9 = J2(h°H)u Jl + fdH + dg : f,g€®(\)[x,y], fi G Q(X)[h], (3.4) 

with deg^j, /, deg^j, g < d and deg h fi < -^ti ■ Moreover, coefficients of f, g and fi 
are of degree at most v + 0(d 3 ) in A and their sizes are bounded by sd , ( ' . 

Proof. It is well known that for any fixed sufficiently generic A one can write de- 
composition (|3.4p with this bounds on degrees in x, y and h and some numerical 
coefficients, see e.g. |G98j . To understand dependence on A, consider Q3.4p as a 
system of linear equations on the coefficients of /, g and fi. Assume that d > n. 
The number of equations (i.e. of coefficients of 9) is (d + l)(d + 2) = 0(d 2 ). The 
number of unknowns is (of coefficients of fi,f and g) is -^V + 0{d 2 ) = 0(d 2 ). 
Coefficients of the left hand side of the equations are polynomials in A, of degrees 
and sizes being 0(d) (coming from fP) and n w correspondingly. By assump- 
tion, on the right hand side are polynomials of degree at most v\ , divided by some 
common polynomial of degree z^, v\ + V2 = v. Their sizes are at most s. Applying 
Cramer rule, we conclude that the coefficients of /, g and /, can be chosen to by 
polynomials of degree v\ + 0(d 3 ) divided by the same common denominator, so of 
degree v + 0(d 3 ), and of sizes sd°( d '. Since the denominators are the same, the 
same bounds hold for the degrees and sizes of /, g and fi . □ 

3.3. Chen homomorphism. Here we prove an analogue of the first claim of The- 
orem[5]for iterated integrals. 

Let U be semigroup algebra corresponding to a semigroup freely generated by for- 
mal variables X\, . . . , X^. We denote the units of the semigroup and of 7ri(r#,p ) 
by e. Let also U be a completion of U in Krull topology corresponding to the 
maximal ideal m = (X%, . . . , Xjv). 

Definition 11. Define the Chen homomorphism ip : tti(Th,Po) —* U as 

tp(S)=e+ ^2 f Wtj. •••Ui k X il ■■■X ih , (3.5) 

where summation is over the set of all non-empty words in alphabet w;. 
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Let j : U — > U/m K+1 U be the natural homomorphism. Define epic as the 
composition ifK — j K f : ^i(Th,Po) — > U/m K+1 U. 

One can easily show that tp (and therefore ipx) is a group homomorphism to the 
set of invertible elements of U (of U /m K+1 resp.), see [H]. 

Note that the space U/m K+1 U is finite-dimensional, and has standard basis of 

monomials {X ix X i2 X ik ,0 < k < K,l <ij < N}. We claim that the image of 

tpx spans U/m K+1 U. 

Lemma 1. Let A— be the set of products of length at most K of the generators 
Sj ofiri(T H ,Po)- The set {j K <p(S),S G A^ k } is a basis ofU/m K+1 U. 

Proof. Claim of the proposition holds simultaneously for all base in H 1 (Th,Z), so 
we can assume that {[Si]} C Hi(T H7 Z) is a dual basis to {[uj t ]} C i/ 1 (r ff ,Z). We 
have tp(e) = e. This implies the statement for K = 0. 

If K = 1, then <p(Si) — ip(e) = X, t + m 2 U, so X { is in the span of j 1 ip(A~ 1 ). 

For K > 1, we see from the previous equality that 

X i± X i2 >>-X ik = (tpiSiJ - tp{e)) (<p(6i h ) - <p(e)) + m K+1 U, 

and, since tp is homomorphism, the right hand side is a linear combination of 
elements of {(p(S),S G A- fe } (mod m K+1 U). 

So {(Pk(S),S G A- fc } spans U/m K+1 U, and, by cardinality reason (here we use 
that 7Ti(r#,po) is a f ree group), is a basis. □ 

3.4. Construction of the horizontal section. Abelian integrals are iterated 
integrals of length 1. The direct analogue of the matrix Si of (|3.1|) for iterated 
integrals of length at most K is the matrix 

S K {H) = \<f Ui,...^, j r ,i s = l,...,N\ (3.6) 

of iterated integrals of length at most K of basic forms uij G B over the cycles 
S = 5j 1 ...Sj k G A- (we adopt convention L = 1). We call these integrals the 
basic iterated integrals. 

The iterated integrals depend on the choice of the base point of tti (Th , Po) ; so we 
choose po as one of the points of intersection of {H = 0} with the line a = {x — 0} 
(generically, there are n+1 such points). Columns of are Sk are just the coordinates 
of (Pk(S), S G A- a ' written in the standard basis of U/m K+1 U. 

For a generic H for all H sufficiently close to H the pairs \T^,po{H)\ are 

diffeomorphic to (Th,Po(H)) by a diffeomorphism close to identity. This diffeo- 
morphism is unique up to isotopy, so we can identify the fundamental groups 
7Ti (T jj,po(H)) . This means that any path S G tti(Th,Po(H)) can be continu- 
ously extended to a family S(H)defined in some neighborhood of H . Therefore 
Sk can be extended holomorphically to some neighborhood of H, and, by analytic 
continuation, to a multivalued matrix function holomorphic on some Zarisky open 
subset of PC n+ i[x, y]. 

Lemma [T] claims that Sk non-degenerate for a generic choice of H. Therefore 
near generic H the matrix Sk describes a basis of sections of the trivial vector 
bundle P<C n+1 [x,y] x U/m K+1 . 

Our goal is to explicitly write coefficients of the connection for which the ma- 
trix Sk is a basis of horizontal sections. We construct this connection locally in a 
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neighborhood of some generic point H £ PC„+i[x, y]. The coefficients of the con- 
nection matrix J7^- = dSx S^- turn out to depend rationally on H and the point 
of intersection p £ {H — 0} n a, which is an algebraic function of H. 

To eliminate the algebraic multivaluedness of Qk we lift the bundle and the 
connection to the corresponding algebraic cover. Namely, for any n > we define 
B n to be the product 

B n = PC n+1 [x,y]xCP 1 (3.7) 

of the space of all polynomials of degree at most n + 1 vanishing at (0,0), and of the 
line a = {x = 0}. Define mapping ev : B n — > PC„ + i[x, y] by ev(H, y) = H — H(y). 
Lifting Clx = ev * Qk defines a meromorphic connection on B n x U/m k+1 . We prove 
that the resulting connections matrix is rational on B n and satisfies the conditions 
of Theorem [U 

3.5. Differentiation of iterated integrals. Our main tool in construction of the 
connection is a formula of differentiating of integrals, a version of the Gelfand-Leray 
formula for non-closed paths. We follow closely |G05] , 

Let R be a functions holomorphic in some open set W C C 2 , and assume that 
its non-critical level {R = 0} is smooth and intersects transversally the line a — 
{x = 0} at point po(0). Choose a path 5 lying on {R = 0} and starting from p$(0) 
with endpoint Pi(0). For any point p in a neighborhood of Pi(0) we can define a 
path S(p) close to (5, lying on {R = R{p)} and joining p and the point po(p) of 
transversal intersection of {R = R(p)} H {x = 0}. 

Proposition 3.3. ( \G05[ Lemma 2.2]) Letui be a differential 1-form holomorphic 
in W . Then for the integral L, > uj, the following equation holds 



d bJ = / -T5)dR + uj-(c7oR)*uj, (3.8) 

J Sip) \Js(p) 

where a : (C, 0) — > {x = 0} is the parameterization of {x = 0} by values of R: 
a(t) = {P = i}n{x = 0}. 

Assume now that the initial path 6 — 6(0) is closed, and the endpoint p of the 
path 6(p) varies on a. As in Definition [U denote the resulting nest of cycles by 6(t), 
where t — R(p). Let u>\, ...,u„ be differential 1-forms holomorphic near 8q. Assume 
in addition that the pullbacks (a o R)*Ui — for the transversal line a. 

Proposition 3.4. The following equation holds: 






dujj 



-p... W „= L ^l... W i-l^ «*! 



Wl...Wi_l — - — u> i+2 ...u> n (3.9) 

art 



i=\ 



Proof. Let us denote r)i = Wi . . . u>i and 9j = u)j . . .u) n . Denote 

(pi(p) = UJi...U) n - 

Jpo 

f7l+l = 1 



ON THE NUMBER OF ZEROS OF MELNIKOV FUNCTIONS 



Also let us define tpi(q) = j q p%, where 



dR 

Then we have 

f p d(wi(pi+i) f p dtpi+iAUi + ifi+idui . 

Vi-iPi= Vi-i T5 = / *?*-! T5 l<Kn 

f P dbJ n 

Vn-lpn = / Vn-l-TB- 

p Jp ali 



and, by I 

ckpi+i = ipi+idR + uj i+1 ip i+2 (3.10) 

hence 

dfi+i Aw, . u>i A Wi+i 

— 55— = w^-i - — 55— ma 

and then 

»7i-ip. = / Wi+i - / Vi-i r^ — Vi+2 + / rji-x—Bi+i 1 < i < n 

Pa Jpo Jpo "'H Jpo "'ft 



f P f P duj n 

/ Vn-lPn = / Vn-l-T^- 

/ mi>i+i = / ' 



'Po •'Po 

Observe that 



iliPi+i 
So we obtain 

Pi = 2^ Vi^—B^ - 2- / ^-1 dff g '+ 2 
Po j = i •'Po "-"" i=1 •'Po "-"- 

Now assume that p — po, so <$(£) are cycles. We will use Gelfand-Leray formula to 
obtain 

d f f d{uj\ip2 

dtf Ul= 

Hence 



xfvi-Zfum-f-M—VPi 



jJe 1 = Y i jra-1%9^ £ <fvi-i^^6i + 2 



Corollary 3.1. Let us assume that Ui = x^ i y li dx, then 



□ 



— j> 0)1 ■ ■ ■ U> n = ^ j> UJl ...Ui-l-p£u>i + l ...W„ (3-11) 

Proof. True since (r o R)*L0i — and Wj A ujj =0. D 
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3.6. Exact forms in iterated integrals. Let wi,...,w„bea holomorphic differ- 
ential 1-forms, g be a holomorphic function in a domain V and 8 C V be a path 
connecting points po and p\ . Our goal is to express iterated integrals involving the 
exact form dg in terms of iterated integrals of smaller length. 

Clearly 

pi 

d 9 = g(pi) -g(po) 



For iterated integrals of length greater than 1, integrating by parts and using (|3.10[) 
gives 



i rPi 

{dg)ui ...cj n = / (dg)(fi = g<pi 

'Po Jpo 

But dR = on level curve, so we have 

rPi 
(dg)u>i ...u n = g(pi) / wi...w„ 
•/po 



gtyxdR + ui^). 



in 



(gui)w2 ■ ■ • w„ 



/'o 



Next, 



/■pi / r? /-q 

Vi(dg)0i+i= / ^(3(9) / ^+1- / (.9^+1)^+2 

Po •'Po \ •'Po •'Po 



Hence 

fin 



/'o 



/■Pi rPi 

Wi . . . Wi(<ig)u; l+ i . . . w„ = / wi . . . (w^)w i+ i . . . w n - / wi . . . Wi(ffWi + i) ...u„ 

Jp Jpo 

(3.12) 



And the third formula: 
Pi rPi 

wi...u n (dg)= / wi...w n (sr(g)-fl'(po)) 



pi /-pi 

wi . . . (u n g)-g(po) / wi 

Po Jpo 



3.7. Construction of Picard-Fuchs system. Let H be a Hamiltonian of degree 
n + 1, which we can write in multi- index form 



A = (Ay) £ PC„ + i[x,y]. 



if = ^ AtfsV 

0<i+j<n+l 

We assume that £> = {a;' y-'da;} = {cj;,Z = l,...,n 2 } form a basis of the Petrov 
module Pjj, and that the curve {H = 0} is smooth and intersects the line a — {x — 
0} transversally. We will compute the connection matrix £Ik locally near H, and 
then, by analytic continuation, this expression will be valid everywhere. 

Let 8 C {H = 0} be a cycle with an initial point at p(H) € Sda, and consider the 

vector of coefficients of <p(8) of (J3.5I) in the standard basis {X^ X.- Lk , k = 1, ..., K} 

oiU/m K+1 : 

( 1 \ 



7(A) 



h 
h 

\In k ) 



where A = {A Q }| Q |< n+1 . We assume that the integrals Ji, . . . , I^ K are ordered by 

length, i.e. Ij = J r\\ . . . % if and only if N^-i < j < Nk, where Nk — dim U/m k+1 . 

Our first goal is to provide an analogue of Proposition 13 . 21 for iterated integrals. 



m . 
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Proposition 3.5. If rji 6 B, . . . ,r]K € B, and 9 is a differential 1-form of degree 
d, then, for any 1 < % < K + 1, 

Arii... rn-i9r)i ...n K = ^ hjlj, h 3 e C(A)[p], (3.13) 

•^ j=l 

where p — 5(0) is the initial point of the cycle S. Degrees of hj in p do not exceed 
d + 0{nK). 

Moreover, if 9 is defined over Q(A), its degree and size do not exceed v and 
s correspondingly and d > n, then the the degrees and sizes of hj do not exceed 
v + 0(K 5 d A ) and s(Kd)°^ d ^ respectively. 

Proof. Using Proposition 13.21 we can write 

i ?7i . . . %_!% ...n K = j> ?yi . . . %_! ^{fj ° E)uj + dg J r)i . . . r} K 

V'- 1 / (3.14) 

rji ...Tfi-iUjTh ..-m + j> m ■■■rii-i(dg)rii ■■■Vk, 

3=1 

since 6 C {H = 0}. The latter term can be rewritten as a sum of two iterated 
integrals of length K, using the equations of N3.61 or, if i — 1 or i = K + 1, as a 
sum of an iterated integral of length K and g(p) § s rji . . . t]k, a multiple of a basic 
integral. This proves (by induction) the formula (J3.13I) . 

Bounds on the size and degree follow from this inductive proof and bounds of 
Proposition 13.21 Note that iterated integrals obtained on different steps of the 
induction have different lengths. First, one can see that deg p hj < d + J2 deg rn < 
d + 2nK, as dependence on p appears only in the remainder g(p) above and their 
degrees are bounded by d + 2nK. 

One can show that all denominators of the polynomials fi,g obtained on al 

inductive steps are factors of A = I Jlfci" A? ) ' w here Aj is the determinant of 
the system Q3.4p corresponding to a form 9 od degree j. We get deg A = 0(K 5 d A ) 
by Proposition 13.21 Therefore, multiplying 9 by A, we can assume that all fi,g on 
all steps of induction are polynomial in A. Denote by b(K + l,d,v) the degrees of 
their coefficients. Then 

b(K + 1, d, v) < maxjV + 0(d 3 ), b(K, d + 2n,v + 0(d 3 ))} 

(polynomiality allows to replace sums of degrees in (|2.3[) by maximum of the de- 
grees). This implies b(K + 1, d, v) < v + 0(K 4 d 3 ). Adding the degree of A, we get 
the required estimate. 

Similarly, assume that the polynomials fi,g obtained on all inductive steps are 
polynomials, and denote by s(K + 1, d, s) the size of the polynomials hj. We have 

s(K + l,d,s)< max |sd° (d3) , 2s(K, d + 2n, sd°^)\ , 

where factor 2 appears due to (|3.12[) . This implies that 

s(K + l,d,s) <s{Kd)° {Kid3) . 

Now, size of A is (Kn)°^ K ' n \ so applying the previous estimate to A9 (and 
remembering d > n), we get the result. □ 
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We can prove more general statement: 

Proposition 3.6. Let 9x,...,6k G A 1 (C 2 ) ® C(A) be 1-forms of degree at most d, 
and of degree in A at most v. Then 

N K 

i...0 K = J2 h ^^ hj€C(X)\p], (3.15) 

with degrees of hj in X,p bounded by vd ^ >. 

Proof. Indeed, decomposing Oi = XLee f iuj °Hlo + fidH + dgi as in Proposition ^. 2| 
we see that 

h...0K = J2\f[ fa® ) j> <I>0)"4{K) + .», (3-16) 




vi=i ' Js 



where summation is over all mappings <f> : {1,...,K} — > B and the dots denote 
(n 2 + 1) K — n 2K iterated integrals of length K with at least one exact form dgi. 
These can be represented as iterated integrals of lesser length, and the result follow 
by induction. 

To estimate the degrees and sizes of the coefficients in (J3.13I) . note that the 
degrees in A of Y[ i=1 jWi) are Kv + 0(Kd 3 ) by Proposition 13.21 The remain- 
ing (n 2 + 1) — n 2K terms can be rewritten by formulae of J j3.6[ as at most 
2 ((n 2 + 1) K — n 2K ^j iterated integrals of length K — 1 with coefficients being ra- 
tional in A and polynomial in p, of degree at most Kv + 0(Kd 3 ) in A and at most 
d in p. Under the integral sign stand some tuple of basic forms, dgi-s and product 
of some gi with one of them. So these are forms of degrees at most 3d. Since the 
degree of the sum of rational functions is at most double their total sum, we get 

b(K, d,v)<A {Kv + 0(Kd 3 ) +b(K-l, 3d, 2v + 0(d 3 ))) ({n 2 + l) K - n 2K ) , 

where b(K , d, v) denote an upper bound for the degree. This implies that b(K, d, v) < 
vd°^\ □ 

However, in order to prove Theorem |4j we will need only dependence on Aoo = —t 
and on p, which are much better: 

Proposition 3.7. Assume that the forms 9i are independent on Aoo- Then the 
coefficients hj in (|3.15j) are polynomial in p, Aoo > an d f or their degrees in p, Aoo we 
get _ 

degif deg Aoo hj + deg p hj + ^ deguj jk = ^ deg0*. (3.17) 

Proof. Polynomiality follows from Proposition 13.21 Therefore the bound can be 
computed by asymptotics at infinity, as t — > oo. For homogeneous forms and 
generic homogeneous Hamiltonians counting homogeneity degrees we get 

deg H deg Aoo hj + dcg p hj + ^ deg U3 ih = ^ deg U 

where Ij = J Wj 1 ...uij m . In general case, this becomes an inequality. D 

Proposition 3.8. Let, as before, H be a Morse-plus polynomial such that the curve 
Th = {H = 0} is smooth and intersects transversally the line a = {x — 0}. Let 
SrW be the Nk X Nx-matrix valued function defined in a neighborhood of H as 
in (|3.6|) . Then 

dS K (X) = n K (\,p)S K (\), (3.18) 
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with Nk x Nk matrix o/f2.R-(A,p(A)) of one-forms on PC n +i[x,y] with coefficients 
being rational functions of A and polynomial in p(X). Here p = po(X) is a starting 
point of integration, p G Th D a . 

The coefficients of Qk{X,p) have degree in p at most O(nK), and their degrees 
in A and sizes are at most 0(K 6 n 7 ) and (Kn) ( K n '. 

Proof. According to Proposition 13. II for any tuple {»?j}|£i, with rji G B, 

^fm...VK = -^Jm...TH-i(j£)fH+i.--VK. (3-19) 

Our goal is to express the Gelfand-Leray derivatives j^- as rational combinations 

>rni Hduji. Differentials 
basis of A 2 (C 2 )/dH A A^C 2 ), so one can write 



of forms ujj. Consider the form Hduji. Differentials of the elements of B form a 



Hduji = y aijduij + rji A dH, i = 1, ..., n . (3.20) 

i=i 

Multiplying by the monomial §t~, dividing by dH and decomposing -w^-r/i as in 
Proposition 13. 21 we get, modulo some multiple of dH, 

H ^=J2 a ^ + [j2 b U H )^+ftdH + d~gA , z,j = l,...,n 2 . (3.21) 

(we follow closely [BNY10, A.3]). Multiplying both sides of this system of equations 
by inverse A of the matrix A = {ay} we get 

where coefficients qfAH) are polynomial in iJ with coefficients being rational func- 
tions of A. 

Integrating over 5 C {H — 0}, we get 

a /> K n 

O 



-q^- pi-lJf = -^^ 9y(0) $??!•■• %-i(Wj)»?i+i 



«=i i=i 



2^ / m---Vi-i(dgThi+i---Vn 



(3.23) 



Using Proposition 13.51 we can express the integrals from the right hand side of the 
latter equation as a combination of basic ones. 

Now, let us estimate the degrees and sizes of the coefficients of £Ik(X,p). The 
degree in x,y of the form rji is equal to the degree of Wj, i.e. is at most In. To 
find the degrees in A and sizes in (|3.20p note that this is a system of 0(n 2 ) linear 
equations on a%j and coefficients of rji of degree 1 in A and of size O(n). Therefore 
dij and rji can be chosen of degree 0{n 2 ) in A and of size n oijl \ 

Bounds of Proposition 13.21 imply that the degrees deg H bfj < 3 = 0(1), degree 
in x, y of gf is bounded by 0(n), and degrees in A and sizes of bf -, gf are bounded 

above by 0(n 3 ) and n°( n ) respectively. This implies that degrees and sizes of 
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&ij are bounded by 0(n 4 ) and n°( n > respectively, so qfj are polynomials in H 

of degree 0(1), of degree 0(n 6 ) in A and of size n°( n '. The functions gf are 
polynomial in x, y of degree 0{n), and their degrees in A and sizes are 0(n 7 ) and 
n o(n ) correspondingly. Bounds of Proposition 1331 now imply that the coefficients 
of f2i<-(A,p) are rational functions of p, A, polynomial in p of degree at most 0{nK) 1 
and of degrees in A and sizes at most 0(K e n 7 ) and (Kn)°( K n >. □ 

3.8. Changing the variables (lifting). The coefficients of the connection (|3.18p 
depend algebraically on A G PC n +i[x, y], since the base point p = p(X) of the 
fundamental group is not defined uniquely. Let ev : B n — > PC n +i[x,y] be the map 
ev(H,y) = H — H(y), where B n was defined in (|3.7[) . Let Sk = Sk ° ev be the 
lifting of the matrix Sk to B n . The coefficients of the pulled-back connection 

dS K = &kS k , &k = ev* n K , (3.24) 

on B n x U/m K+1 are rational one- forms on B n . 

Proposition 3.9. Degree and the size of the matrix Q-k are bounded as 0(K 6 n s ) 
and (Kn) ( K n ' correspondingly. 

Proof. Degree of the mapping ev is n + 1, and it has coefficients equal to or 
1. Therefore degree of flu is at most 0(K 6 n 8 ). Sizes of coefficients of £Ik will 
not exceed sizes of coefficients of $7#-(A,p) multiplied by the size of ev raised to 
deg A fl K (\p), i.e. (K n )0(K 8 » 5 ) n 0(KV) = (^ n )0(KV)_ D 

4. Properties of the system 

4.1. Quasi-Unipotency. 

Proposition 4.1. Connection (|3.24[) is quasi-unipotent and regular. 

Proof. From p.23[) we see that derivatives of an iterated integral are linear com- 
bination of iterated integrals of smaller or equal length. This means that £Ik is a 
lower-block-triagonal matrix: 

/O \ 

0ii 

* 022 

* ••• 

where each block 0^ is ki x ki matrix corresponding to the integrals of length 
exactly i Note that ©n is just the pull-back by ev* of the matrix of the Gauss-Manin 
connection for Abelian integrals, so has quasiunipotent monodromy by [K81 . 

Since VLk is lower-block-triagonal, it preserves the flag T = {0} = J^k+i C 
Tk C ... C Tq = U/m K+1 , where T\ = m i /m K+1 , so any monodromy operator 
corresponding to £Ik preserves this flag as well, i.e. is lower-triangular with blocks 
Ma on diagonal. Therefore its eigenvalues are just the eigenvalues of the mon- 
odromy operators Mu corresponding to &u, the connection induced by Qk on the 
factor-bundle with fiber J^/J^+i = m t /m 4+1 . 

Lemma 2. For all 1 < k < n, Qkk = P(©j = i 0n)P _1 f or some permutation 
matrix P , where we use the notation of Kronckcr sum: A® A = A® I -\- 1 ® A. 
Correspondingly, M kk = PMf*P~ x . 



n K = 
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The first claim is just a way to say that the differentiation of iterated integrals 
satisfies Leibnitz rule, up to iterated integrals of lesser length. This can be seen 
from p. 191) . Now, derivation of the products of Abelian integrals satisfy the same 
rule, so horizontal sections of Qkk are described by these products, and monodromy 
operators of are just tensor powers of monodromy M\\ of Abelian integrals. This 
proves quasiunipotency of VLk since tensor powers of quasiunipotent operators are 
quasiunipotent. 

Regularity of (Ik follows from regularity of Q k k and the fact that semidirect 
product of regular connections is regular, see [D] , 

□ 

5. Proof of Theorem [2] 

Let 8 C {H = 0} be a cycle and U its small neighborhood, and let M^ be 
the first non-zero Melnikov function defined as in (jl.31) . Recall the construction 
expressing Mr- as a polynomial in iterated integrals of the perturbation form u> 
and its Gelfand-Leray derivatives up to order K. 

Definition 12. A real analytic 1-form a G A 1 (C/) is relatively exact with respect 
to the integrable foliation T = {dH — 0} in a domain U, if 

a = h-dH + dg, h,ge 0{U) (5.1) 

Clearly, the integral of a relatively exact form a along any closed oval on any 
level curve {H — z} C U, vanishes: 

V oval 6 C {/ = z} I a = (5.2) 

Define the sequence of real analytic 1-forms u\, k>2) • • • > Wfe as follows: 

(1) (Base of induction). u>\ = lo is the perturbation form from (|1.2[) 

(2) (Induction step) . If the forms u)\, . . . , u)j are already constructed and turned 
out to be relatively exact, then uij — hj ■ dH + dgj. In this case we define 

uij+i = —hjio (5-3) 

Theorem 6. ]IY\ Theorem 26.7] If Wk, k > 2, is the first not relatively exact 
1-form in the sequence W\, . . . ,Wk-ii<^k, constructed inductively by (I5.3|) . then 



M k (z) = -<h co k (5.4) 

J{H=z} 



Evidently, the functions hj can be restored as hj = — J -iri- , 



d ( f d ( ( f duj 



1 i dH \ J dH V V J dH 

Denote by <f> the algebraic function H\~ x o H of x, y which maps the point [x, y) 
to the (one of d + 1) point of the transversal a — {x = 0} lying on the same level 
curve of H as (x, y). In other words, p = 4>(x, y, A). 

Lemma 3. hj in (I5.3[) is a linear combination of iterated integrals of differential 
one-forms with coefficients polynomial in x, y. The coefficients of this combination 
are rational in X,p. 
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Proof. We have ho — 1, which is of the required form trivially. We proceed by 
induction on j. Assume that hj is a finite sum of terms of the type R(X,p) J 9\...6k, 
where Oi are polynomial in x, y one-forms, and R is a rational function. Now, 
hj+x = J ,j, , and it is look what happens with just one term of this sum. 
Differentiation, according to 13.31 gives 

^('fl(A,p)«y» 1 ...«^ = |2(ifi;(rt)" i y^i..A+ 

nJp,..h + R±J.»,..^..A + (5 . 5) 

i— 1 



fc- 



' A ^+i a , d f /> Ok^ A (o o H) * 6 k 



r ^o l .. u -^^...e k + R 



dH J dH 

The first term is clearly of the required type. 

2 

Denote by m(H) the product Yi7=i (H ~ c i)i wnere c i are critical values of H 
(repeated if multiple). It is well known, see |G05] [Prop. 2.4]. that m(H) lies in 
the Jacobian ideal of H, so the operator m j H preserves polynomial one- forms. 
Therefore 



^ 6 0_iMi±L = R fj^ a - m(F)^A^ +1 



^-^ dH m(H(p)) J ^ dH 



R sh f a rn(H)d6i 

k = ,tti w / . I wy l- 



(5.6) 



m(H(p))^J " dH - 

all terms in (|5.5p except the last one are of the required type. 

Finally, [a o H) *9 = cf)k(p)dH, where 4>k{p) = dH()(d) 1S a ra tional function of 
p, so the last term of ()5.5[) is also of the required type. □ 



Lemma 4. For a form uj of degree d > n in x,y, we have 

N K 

M K =J2hiIi, 

i=i 

where hj depend rationally on p and has degree at most 2 ( K 'dn 6 in p. 

Proof. In the inductive step (|5.5j) one iterated integral of forms of degrees at most 
dk with coefficient R of degree Vk generated 0(k) iterated integrals of forms of 
degrees at most 2dk + 0(n 3 ). The coefficients of these new integrals are obtained 
from R by a combination of a differentiation and division either by 7?|^(p) or by 
m(H(p)) or just by division by one of these polynomials. Applying these operation 
K times we can increase degree of R by at most 0(Kn 3 ). Summing together, we get 
a representation of Mk as a sum of 2°( K ' iterated integrals of forms of degree at 
most 2°( K } (d + 0(n 3 )), with rational in p coefficients of degree at most 0(Kn 3 ), 
with common denominator of degree 0(Kn 3 ). 

Applying Proposition 13.71 we represent each of these iterated integrals as in 
(|3.15[) . and the coefficients of these representations have degrees in p at most 
2°(- f 0(d 4- n 3 y Summing these representations together, we arrive to the state- 
ment of the Proposition. □ 



ON THE NUMBER OF ZEROS OF MELNIKOV FUNCTIONS 



17 



Multiplying by the common denominator, we see that the estimate of Theorem [2] 
follow from the following 

Lemma 5. Linear combination X)i=i hi^i °f iterated integrals of length at most 
K , with coefficients hi being polynomial in p of degree at most /i, has at most 
exp (exp y/j, ^'n W)J zeros on each line A x C p . 

Proof. We can construct in a standard way a connection whose horizontal sections 
are described by these functions. More exact, the connection will be a Kronecker 
sum of Qk with the connection whose sections are polynomials in p of degree at 
most jjl. Therefore it will be of dimension (J,Nk = fJ.n°( K \ of degree 0(K 6 n 8 ) and 
of size /i(Kn) ( K n '. Theorem 0] then gives the required upper bound. □ 

Substitution of /j = 2°^ K ^dn e from Lemma |4] provides the estimate of the The- 
orem [2] 
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